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Abstract. In 2003 a previously unpulsed Einstein and ROSAT source cataloged as soft and dim (L^ of few xl0 33 ergs _1 ) 
thermal emitting object, namely XTE J1810-197, was identified as the first unambiguous transient Anomalous X-ray Pulsar. 
Two years later this source was also found to be a bright highly polarized transient radio pulsar, a unique property among both 
AXPs and radio pulsars. In September 2006 Swift Burst Alert Telescope (BAT) detected an intense burst from the candidate 
AXP CXOU J164710. 2— 455216, which entered in an outburst state reaching a peak emission of at least a factor of 300 
higher than quiescence. Here, we briefly outline the recent results concerning the outburst phenomena observed in these 
two AXPs. In particular, XTE J1810-197 has probed to be a unique laboratory to monitor the timing and spectral properties 
of a cooling/fading AXP, while new important information have been inferred from X-ray and radio band simultaneous 
observations. CXOU J164710.2— 455216 has been monitored in X-rays and radio bands since the very beginning of its 
outbursting state allowing us to cover the first phases of the outburst and to study the timing and spectral behavior during 
the first months. 
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INTRODUCTION 

At the beginning of the X-ray astronomy era the study of 
the X-ray variability (both of the flux and timing proper- 
ties) of Cen X-3 allowed astronomers to unambiguously 
assess the binary nature of the source and to identify the 
accretion of mass, flowing from the companion onto a 
rotating neutron star, as the main mechanism to produce 
X-rays 01 01- Since then, several classes of high energy 
sources hosting a compact object have been identified 
and their variability, observed on different timescales, 
used to test models and/or to study the emission mech- 
anisms as a function of flux (while keeping fixed other 
parameters, such as the distance and the geometry of the 
system). For the class of accreting neutron stars this ap- 
proach was successfull in confirming, among others, the 
presence of a centrifugal barrier to accretion, testing the 
dependence of the period derivative with respect to the 
source luminosity, and studying the change in the pho- 
ton propagation direction as a function of the accretion 
rate by means of the pulse shape changes [fan and pen- 
cil beam model; as an example see B |4j, |5J . Isolated 
neutron stars are relatively constant sources, making the 
above approach unreliable. However, there are two small 
classes of isolated neutron stars which show spectacular 
events during which their luminosity may change up to 
10 orders of magnitude on timescales down to few mil- 



liseconds. These objects are better known as Anomalous 
X-ray Pulsars (AXPs; 10 objects plus 1 candidate) and 
Soft y-ray Repeaters [SGRs; 4 objects plus 2 candidates; 
for a review see|6J]. It is believed that AXPs and SGRs are 
linked at some level, owing to their similar timing prop- 
erties (spin periods in the 2-12 s range and period deriva- 
tives P in the 10 13 — 10 11 ss _1 range). Both classes 
have been proposed to host neutron stars whose emis- 
sion is powered by the decay of their extremely strong 
magnetic fields [> 10 15 gMM- 

Different types of X-ray flux variability have been dis- 
played by AXPs. From slow and moderate flux changes 
(up to a factor of few) on timescales of years (virtu- 
ally all the object of the class), to moderate-intense out- 
bursts (flux variations of a factor up to 10) lasting for 1-3 
years (IE 2259+586, and IE 1048.1-5937), to dramatic 
and intense SGR-like burst activity (fluence of 10 36 — 
10 37 ergs) on sub-second timescales [4U 0142+614, 
XTEJ1810-197, 1E2259+586 and IE 1048.1-5937; 
see |9j, for a review on the X-ray variability]. The first no- 
table recorded case of flux variability was the 2002 burst- 
ing/outbursting event detected from IE 2259+586, the 
only known event in which a factor of ~ 1 persistent flux 
enhancement in an AXP was followed (or proceeded) by 
the onset of a bursting activity phase during which the 
source displayed more than 80 short bursts CaEl]. The 
timing and spectral properties of the sources changed sig- 



nificantly and recovered the pre-bursting activity phase 
values within few days, likely due to the relatively high 
luminosity DC level (~ 10 35 ergs~'). However, it was 
only in 2003 that the first transient AXP was discov- 
ered, namely XTE J1810-197, which displayed a fac- 
tor of >100 persistent flux enhancement with respect 
to the unpulsed pre-outburst quiescent luminosity level 
[~ 10 33 ergs _1 ; [II El [II El- Unfortunately, the ini- 
tial phases of the outburst were missed and we do not 
know whether a bursting activity phase, similar to that of 
IE 2259+586, occurred also for this source, though four 
bursts have been detected by RossiXTE between 2003 
September and 2004 April and unambiguously associ- 
ated with XTE Jl 8 10- 197 d]. In 2006 the candidate 
AXP CXOU J 1647 10.2-4552 16 displayed a bursting- 
outbursting behavior with a maximum flux variability of 
> 300, followed by extreme and dail y chang es both in 
the spectral and timing properties lll7lll8[|l9ll . 

These two sources currently represent our best oppor- 
tunity in order to study the evolution of the main spectral 
and timing parameters as a function of flux by keeping 
fixed other parameters (such as distance and geometrical 
angles) otherwise difficult to infer (similarly to the pio- 
neering studies on accreting X-ray pulsars). In the fol- 
lowing pages we will briefly outline (i) the recent results 
obtained from the analysis of the latest 4 years XMM- 
Newton monitoring observations of XTE J18 10-197 as 
it approached to quiescence, (ii) the comparison of X- 
ray and radio emission from XTE J18 10-197 by means 
of two ^8hr-long simultaneous campaigns with XMM- 
Newton and Parkes, (iii) and the results of the first 6 
months monitoring of CXOU J 1 647 1 0.2 -4552 1 6 during 
the first phases of its current outburst by means of Swift, 
Chandra and XMM-Newton in the X-rays and Parkes in 
the radio band. 



XTE J1810-197: FROM OUTBURST TO 
QUIESCENCE 

Since the very first XMM-Newton 2003 obser- 
vations of XTE J18 10-197, carried out approxi- 
mately one year after the onset of the outburst, 
it was evident fl3ll that the source spectral shape 
(two blackbodies with kT=0.29±0.03keV and 
R BB w5.5km, and kT=0.70±0.02keV and R BB ssl.5 km; 
Lx ^5xl0 34 ergs _1 in the 0.5-10keV range) was sig- 
nificantly different from that serendipitously recorded 
by ROSAT in 1992 (one BB with kTwl60eV and 
R BB «10km; extrapolated luminosity in the 0.5-10keV 
range of Lx ^7x 10 32 ergs~' and for a distance of 
3.3 kpc). Moreover, the source showed a 5.54 s pulsation 
with a pulsed fraction of nearly 45% during outburst, 
while an upper limit of 24% was inferred from the 



ROSAT data. The above issues originated a number 
of important questions awaiting for an answer: Is the 
soft BB component detected by XMM-Newton evolving 
into the quiescent BB component seen by ROSAT ? 
Alternatively, is the emission from the whole surface 
always present ? What happens to the higher temperature 
BB component as the source approaches to quiescence 
? Which is the pulsed fraction level of the source in 
quiescence (if detectable) ? Does the outburst changed 
permanently the timing/spectral properties (such as the 
pulsed fraction, the flux and temperature or size of the 
quiescent BB component) of the source ? 

In order to try answering to the above questions we re- 
duced all the archival (6) and still proprietary (2) XMM- 
Newton observations [for the details see |20j 12111 and fit- 
ted the eight spectra all together. All the spectra have 
been rebinned in order to ensure that each background- 
subtracted spectral channel has at least 25 counts, and 
that the instrumental energy resolution is not oversam- 
pled by a factor larger than 3 [22, ; indeed the correct 
application of the above rules prevents artificially low 
(good) reduced X 2 s.]. In particular, we can outline the 
obtained results as follows: 

2BB model: By extending the spectral recipe outlined 
by E l23ll we applied the two BB spectral fit analy- 
sis to the fading phases of XTE Jl 8 10- 197 until March 
2007 when the flux source was ~1.2 times above the 
pre-outburst level (reduced x 2 ~1.23 for 975 degree of 
freedom, d.o.f.; N//=0.58±0.02 x 10 22 cm" 2 ). While the 
soft BB component smoothly approaches to that in qui- 
escence (see Figure Q] left panel, 2 nd and 3 rd plots), we 
note a number of ambiguities difficult to account for by 
means of simple assumptions. The hard component BB 
radius is not monotone and it increases after 2.5 years of 
smooth decrease (left panel, 4 th plot) while the tempera- 
ture approaches to that of the soft BB in 2003 (left panel, 
5'' 1 plot). Moreover, none of the spectral parameters or 
components is able to account for the flattening, at the 
25% level, showed by the pulsed fraction evolution (left 
panel, V plot) & 

3BB model: The addition of a further BB com- 
ponent gives a better fit (reduced % 2 ~1 . 15 for 973 
d.o.f.; N H =0.70±0.02xl0 22 cm- 2 ; F-test probability 
gives 7.3(7) though not yet satisfactory. Notably, the fit 
gives parameters and flux (for the coldest BB) which 
are virtually equal to those inferred in quiescence. Even 
more interesting, the hottest BB components show a 
nearly constant evolution of the temperatures (see Fig- 
ureQ] right panel, 3 rd ad 5 th plots), leaving the radii 
as the only variable parameters to account for the de- 
caying phases of the outburst (right panel, 2 and 
4 th plots). Since September 2006 the hottest BB is 
not anymore needed to fit the spectra (upper limit of 
~5 x 10 _15 ergcm _2 s _1 ). In this scenario, the already 
mentioned flattening of the pulsed fraction might be eas- 
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FIGURE 1. Evolution of the spectral parameters for the 2BB (left panel) and 3BB (right panel) models together with the pulsed 
fraction as a function of the observed 1-10 keV band flux. Central Panel: schematic view of the assumed scenario with marked the 
emission regions of the hottest BB components. 



ily accounted for by the disappearance of the hot BB. 
A pulse phase spectroscopic analysis do not show any 
phase lag between the two highest temperature BBs ll2lll . 

Further components: During the first two XMM ob- 
servations (2003-2004) the spectral fit residuals clearly 
suggest (at a 3.2c confidence level) the likely presence 
of an additional hard component above 7-8 keV which 
we are not able to characterize due to poor statistics in 
this band. We can only speculate that might be somewhat 
related to the presence of a hard power-law-like compo- 
nent detected in other AXPs JUS by INTEGRAL and 
which extends up to (at least) 200keV [26]; 

X-ray and radio campaigns: The simultaneous ra- 
dio and X-ray observations of XTEJ18 10-197 carried 
out in September 2006 and March 2007 showed that the 
pulse alignment between the two bands is high and stable 
(see Figure|2]i, while the pulse width is relatively small 
(^0.1 in phase) in the radio, though not unusual among 
radio pulsars $2j[ l28ll . This suggests that the X-ray and 
radio emitting ragions are likely different but nearby (or 
superimposed), the X-rays likely coming from a larger 
area. Moreover, during the first campaign large radio flux 
(~50%) and pulse shape variations have been detected 
which do not correlate with any change (at a few per- 
cent level) of the X-ray timing and/or spectral parame- 
ters. This suggests that the X-ray emission likely orig- 
inates deep in the crust and, more in general, the radio 
and X-ray mechanisms appear different. 

We emphasize that the spectral model used above 
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FIGURE 2. XTEJ1810-197 XMM-Newton PN+MOS 0.5- 
lOkeV and Parkes 20cm light curves, folded to the spin pe- 
riod, carried out during the September 2006 campaign. Super- 
imposed to the X-ray folded light curve is the best sinusoidal 
fit [adapted from l27h . 



(three BBs) is a first attempt to infer the evolution of a 
number of physical quantities while making use of well 
assessed and reliable components. The three BB model 
has the advantage of being model-independent, and of 
minimizing the number of variable parameters during 
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FIGURE 3. Left Panel: phases of the Swift XRT, XMM-Newton and Chandra observations of CXOU J164710.2-455216: a large 
and quick decaying component is clearly present. Note that the first XMM-Newton point (at day -5) would be at the reported phase 
only in the hypothesis that the pre- and post-glitch parameters are similar [for more details see fLUl . Right Panel: the 15-150 keV 
Swift BAT spectrum of the 20ms long burst detected from CXOU J164710.2-455216 on 21st September 2006, together with the 
5ms-binned BAT lightcurve in the time interval around the trigger (inset). 



the outburst (only the radii are changing and both de- 
creasing). With respect to other model recently devel- 
oped [see for example 12911 we carried out a cross check 
with the timing properties and rejected all those mod- 
els/components not able to reproduce the pulsed frac- 
tion evolution. All the above inferred information are 
important in the effort we are currently making in de- 
veloping and using more detailed and complex (but nec- 
essarily model-dependent) modelizations. In particular, 
the simultaneous fit of the spectra and energy resolved 
folded light curves might provide a tool to infer the ge- 
ometry of the surface temperature distribution and to in- 
dependently check the goodness of the assumed spectral 
model(s). Finally, we note that the three BB model fit im- 
plies that the BB component accounting for the emission 
from the whole sourface is almost unpulsed (pulsed frac- 
tion of 9%±2%); this prediction can be easily checked 
once the source will return into the quiescent state. 

CXOU J164710.2-455216: FROM 
QUIESCENCE TO OUTBURST 

On 2006 September 21, the candidate AXP CXOU 
J 1647 10.2 -4552 16 emitted an intense (~10 39 ergs _1 ) 
and short (20ms) burst promptly detected by the Swift 
BAT. Together with the burst, large changes in the tim- 
ing and spectral properties of the persistent component 
were detected and seen evolving during the subsequent 
weeks. In particular, the Swift XRT monitoring (plus two 
proprietary XMM-Newton and two archival Chandra ob- 
servations ) during the first six months since the outburst 



allowed to infer the following characteristics: 

The BAT burst: The prompt event recorded by Swift 
BAT has an exponential time decay T of 3.3±1.0ms (lrj 
confidence level) and the spectrum can be fitted with 
both a blackbody with kT of 9.9+Jj keV and a T of 
1.8 ±0.5 (see Figure[3] right panel). In both cases a flu- 
ence of ~ 10 -8 ergcm~ 2 corresponding to a total en- 
ergy of ~ 2 x 10 37 ergs (for a fiducial distance of 5kpc). 
Compared with the properties of the previously detected 
AXP bursts, the current burst has a duration within 1 a 
from the log-normal distribution average value inferred 
for IE 2259+586, while the fluence is significantly (a 
factor of about 50) larger than the mean iflQfl . There is 
only one burst, out of ^80 detected from IE 2259+586, 
with duration and fluence comparable with that of CXOU 
J164710.2-455216, while a total of three bursts (over 
the whole burst duration distribution) have fluence com- 
parable or slightly larger than that of BAT. 

The phase-coherent timing and the glitch: The 
pulse phase evolution is consistent with the occurrence 
of a large glitch (Av/v ~ 10~ 4 ), the largest ever de- 
tected from a neutron stafl We also detected a quadratic 
component in the pulse phases corresponding to a P = 
9.2(4) x 10~ 13 s s _1 and implying a magnetic field 
strength of 10 14 G (see Figure[3] left panel). The first 1- 



The glitch detection was obtained by minimizing the number of 
variable peaks in the pulse profile. A less significant timing solution 
is feasible and requesting only a P component. We consider the latter 
unlike since it would imply high variability for all the peaks; see fl^l 
for details. 



lOkeV Swift XRT spectrum was carried out ~13 hours 
after the burst detection and showed, in addition to a 
kJ ~ 0.65 keV blackbody (R BB ~ 1.5km) , a T ~ 2.3 
power-law component accounting for about 50% of the 
observed flux (alternatively, two blackbodies with kT s = 
0.50 ± 0.05 keV with R BBs = 3.2 ± 0.4km, and kT h = 
l.l±o jkeV with R BBh = 0.5 ±0.1 km; all the uncertain- 
ties are at the 90% level); 

The flux and pulsed fraction evolution: The flux 
decay of CXOU J 1 647 1 0.2 -4552 1 6 is well described by 
the function F « t a , with index a of -0.28±0.05 (similar 
to the case of the 2002 IE 2259+586 burst-active phase). 
Moreover, we found that the PL component decays more 
rapidly (index a of -0.38±0.1 1; 90% uncertainty) than 
the BB flux (index a of -0. 14±0. 10). The pulsed fraction 
of the 10.61 s pulsations was seen to drop from a value 
of ~ 80% (as recorded by an XMM-Newton observation 
few days before the burst) to ~ 10% few hours after the 
BAT event. The spectral and timing analysis clearly show 
that only the blackbody component is responsible for the 
pulsed flux. In particular, the signal fractional rms as a 
function of time is well fitted by a power-law with index 
a of +0.38±0.11, equal (but with opposite sign) to the 
power-law component decay. 

The quiescent properties: Archival Chandra data 
analysis revealed that the modulation in quiescence is 
100% pulsed at energies above ^4keV and consistent 
with the (unusually small-sized) blackbody component 
being occulted by the neutron star as it rotates. 

Radio observations: Since the onset of the outburst, 
CXOU J 1647 10.2-4552 16 has been routinely observed 
in the radio band from Parkes. No pulsed radio emission 
has been detected so far with an upper limit of 0.04mJy 
at 1400MHz (Ml. 

All these results confirmed unambiguously that 
CXOU J164710.2-455216 is a transient and bursting 
AXP, showing an unusually high pulsed fraction level 
in quiescence. In particular, the comparison of the cu- 
mulative properties of CXOU J164710.2-455216 with 
those of other AXPs which showed a similar behavior in 
the latest years confirmed that these outbursting events 
are more common than previously thought.The way 
and timescales with which the post-burst/glitch timing 
and spectral properties will recover to those in quies- 
cence hold a great potentiality of better understanding 
the mechanisms which rule the outbursts of AXPs. In 
particular, the comparison with the XTEJ18 10-197 out- 
burst evolution might help us in understanding whether 
other parameters play an important role in the observed 
properties. 

Finally, we note the the BAT detection of the bursts 
from CXOU J 1 647 10.2— 455216 opens new prospectives 
for detecting further burst from known AXPs and for 
identifying new AXPs/SGRs with the Swift mission. 



ACKNOWLEDGMENTS 

The complete list of the co-authors is: S. Campana, A. 
Corongiu, J. Cummings, M. Dahlem, M. Falanga, P. Es- 
posito, D. Gotz, S. Mereghetti, P. M. Muno, R. Perna, 
R. Turolla and S. Zane. This work is partially supported 
at OAR through ASI, MIUR, and INAF grants. We 
acknowledge financial contribution from contract ASI- 
INAF 1/023/05/0. 



REFERENCES 

1. Giacconi, R., Gursky, H., Kellogg, E., Schreier, E., & 
Tananbaum, H., ApJ, 167, L67 (1971) 

2. Schreier, E., Levinson, R., Gursky, H., Kellogg, E., 
Tananbaum, H., & Giacconi, R., ApJ, 172, L79 (1972) 

3. Parmar, A. N., White, N. E., Stella, L., Izzo, C, & Ferri, P., 
ApJ, 338, 359 (1989) 

4. Parmar, A. N, White, N. E., & Stella, L., ApJ, 338, 373 
(1989) 

5. Angelini, L., Stella, L., & Parmar, A. N„ ApJ, 346, 906 
(1989) 

6. Woods, P. M., & Thompson, C, Compact stellar X-ray 
sources, Cambridge University Press, Cambridge, pp. 
547-586, astro-ph/0406133 (2006) 

7. Duncan, R. C, & Thompson, C, ApJ, 392, L9 (1992) 

8. Thompson, C, & Duncan, R. C, MNRAS, 275, 255 (1995) 

9. Kaspi, V., Isolated Neutron Stars: From the Interior to the 
Surface, Astrophysics & Space Science, Springer, Vol. 308, 
Issue 1-4, pp. 1-11, astro-ph/06 10304 (2007) 

10. Gavriil, F. P., Kaspi, V. M., & Woods, P. M., ApJ, 607, 959 
(2004) 

11. Woods, P. M., et al., ApJ, 605, 378 (2004) 

12. Ibrahim, A. I., et al., ApJ, 609, L21 (2004) 

13. Gotthelf, E. V., Halpern, J. P., Buxton, M., & Bailyn, C, 
ApJ, 605, 368 (2004) 

14. Israel, G. L., et al., ApJ, 603, L97 (2004) 

15. Rea, N., et al., A&A, 425, L5 (2004) 

16. Woods, P. M., et al. ApJ, 629, 985 (2005) 

17. Muno, M. P., et al., ApJ, 636, L41 (2006a) 

18. Muno, M. P., et al., MNRAS, 378, L44 (2006b) 

19. Israel, G. L., et al., ApJ, 664, 448 (2007a) 

20. Rea, N., Nichelli, E., Israel, G.L., et al, MNRAS, in press, 
astro-ph/0707.3363 (2007) 

21. Bernardini, F, Israel, G.L., Rea, N, et al., A&A, submitted 
(2007) 

22. Ibarra, A., Matt, G, Guainazzi, M., et al., A&A, 465, 501 
(2007) 

23. Gotthelf, E. V., & Halpern, J. P., ApJ, 632, 1075 (2005) 

24. Kuiper, L., Hermsen, W., & Mendez, M., ApJ, 613, 1 173 
(2004) 

25. Kuiper, L., Hermsen, W., den Hartog, P. R., & Collmar, 
W.,ApJ, 645, 556 (2006) 

26. Gotz, D., et al. The Astronomer's Telegram, 953, (2006) 

27. Israel, G. L., et al., ApJ, submitted (2007b) 

28. Lorimer, D. R., et al., MNRAS, 372, 777 (2006) 

29. Guver, T. et al., ApJ, 667, L73 (2007) 

30. Burgay, M., Rea, N., Israel, G., & Possenti, A., The 
Astronomer's Telegram, 903 (2006) 



